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ABSTRACT 

We present a new XMM-Newton spectrum of the Seyfert 2 nucleus of IC 2560, 
which hosts H 2 maser emission from an inclined Keplerian accretion disk. The 
X-ray spectrum shows soft excess due to multi-temperature ionized plasma, a 
hard continuum and strong emission features, from Mg, Si, S, Ca, Fe and Ni, 
mainly due to fluorescence. It is consistent with reflection of the continuum 
from a mostly neutral medium and obscuration due to a high column density, 
> 10 24 cm -2 . The amplitude of the reflected component may exceed 10% of the 
central unobscured luminosity. This is higher than the reflected fraction, of a 
few percent, observed in other Seyfert 2 sources like NGC 4945. We observe an 
emission line at 6.7 keV, possibly due to FeXXV, undetected in previous Chandra 
observations. The absorption column density associated with this line is less than 
10 23 cm -2 , lower than the obscuration of the central source. We hypothesize that 
this highly ionized Fe line emission originates in warm gas, also responsible for 
a scattered component of continuum emission that may dominate the spectrum 
between 1 and 3 keV. We compare X-ray and maser emission characteristics of 
IC 2560 and other AGN that exhibit water maser emission originating in disk 
structures around central engines. The temperature for the region of the disk 
associated with maser action is consistent with the expected 400-1000K range. 
The dumpiness of disk structures (inferred from the maser distribution) may 
depend on the unobscured luminosities of the central engines. 
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INTRODUCTION 



IC2560 is a nearby spiral 



cz 



2925 km s x ) that is classified Seyfert 2 (jFairall 



19861 ) and known t o host 22 GHz H ? Q m aser emission originating in an accretion disk around 
the central engine (jlshihara et a/.ll200ll ). Here we present analysis of a deep (archival) XMM- 
Newton EPIC observation of this active galact ic nucleus (AGN). P revious observations of 
IC 2560 in the X-ray band using ASCA GIS Jlshihara et all hoOlh and Chandra ACIS-S 
(jlwasawa et a/.ll2002l ; iMadejski et a/.ll2006l ) both pointed to a reflection dominated spectrum 
and high obscuring column density. 



IC 2560 is a relatively weak X-ray source, with an observed 2-10 keV flux of 3.3x10 



-13 



erg s 1 cm 2 that originates predominantly from a central source (IMadejski et al\\200w . 



Interpretation of early GIS and ACIS-S observations was constrained by relatively limited 
sensitivity and energy resolution, in particular above 6 keV, where the bulk of the emission 
due to nuclear activity escapes into our line of sight. The EPIC instrument provides ~5x 
larger collecting area at 6 keV and ~20% finer energy resolution (160 eV) than ACIS-S (chip 
3). Although dominated by the central point source, flux detected with EPIC (1.1 kpc pixel -1 
at 1 keV) will include contribution from extended and off-nuclear emission components. As 
such, we note that the nuclear structure in IC 2560 appears to be relatively complex, as may 
be anticipated for a type-2 object. Hubble Space Telescope ACS/HRC ob servations at 3300 A 
detec t patchy emission with multiple peaks within the central ~ 200 pc (IMunoz Marin et al. 



20071 ). presumably due to dust. Relatedly, population synthesis models fit to optical spectra, 
obtained with the ESO 1.5m (3470-5450 A) fo r the central ~ 300 pc, a re suggestive of 



active star formation during the last 1.5-10 Myr (jCid-Fernandes et all 120041 ) . Ground-based 
spectroscopic study has also detected asymmetry in the [OIII] line profile (a blue excess) 
that may be indi cative of outflow at ~ 10 km s _1 projected on the sky plane in a 300x400 
pc central patch (jSchulz fc Henkel II2003I ). Because the inclination of the galaxy is 63°, and 
the inclination of the central engine (as inferred from maser emission) may be even closer to 
edge-on, the actual flow speed may be larger. 

The observed H 2 maser spectrum exhibits emission close to the systemic velocity as 
well as blue and red-shifted co mplexes, offset, approximately symmetrically, by ± 200 - 
420 km s -1 flBraatz et al\ 120031 ). This pattern is an indicator of emission from a rotating 



edge-on disk structure (e.g., NGC 4258: iNakai et all (Il993l )). Ilshihara et all (|200ll ) reported 
line-of-sight accelerations for the "systemic emission" and mapped its emission distribution 
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using Very Long Baseline Interferometry (VLBI). Based on the assumption that the emis- 
sion traced a thin, rotating, circular disk that is observed edge-on, they used the velocities 
of high- velocity emission to estimate disk radius and dynamical mass. Follow-up VLBI ob- 
servation by Greenhill et al. (2008, in preparation) enabled mapping of both systemic and 
high-velocity emission, thus tracing the disk structure, demonstrating Keplerian rotation of 
material around the central engine in a relatively thin distribution, and enabling more certain 
modeling of the accretion disk geometry and estimation of the central mass. The observed 
molecular gas lies at radii 0.08 - 0.27 pc and the inferred enclosed mass is 2.9xlO 6 M with 
~ 20% accuracy. The thinness of disk (h/r <Cl) as well as presence of H2O emission suggest 
disk material that is relatively cool (<10 3 K). 

Due to generally high extinction toward Seyfert-2 central engines, X-ray and radio ob- 
servations together are especially useful in the study of the physical processes that occur at 
radii <1 pc. Measurement of X-ray spectra enable estimation of intrinsic luminosity, obscur- 
ing column density, temperature and abundance. Time monitoring facilitates estimation of 
size scales for high-energy phenomenon. VLBI observations of H2O maser emission enable 
mapping of accretion disk geometry and orientation to the line of sight, and estim ation of or- 



bital r adii for molecular material, central engine mass and Eddington luminosity (jGreenhill 



( 120071 ) and references therein). The presence of maser emission implies the existence of a 
reservoir of relatively cool material and at the same time, a heating mechanism to maintain 
population inversion. The ne cessary energy may be i mparted by viscous heating within the 



disk (e.g., flDesch et q/.lll998l )) or external irradiation (INeufeld et al\\ 1994 1 Watson fc Wallin 
19941 ). However, in either case, survival of molecular gas demands a shielding column from 
UV and X-ray emission generated by, e.g., the central engine and coronae. This column 
may arise in outflowing material, disk surface layers, or disk material at small radii, thus 
providing additional quantitative constraints on overall structure around the central engine 
and energetics. In section 2, we discuss the data reduction of the EPIC data, while section 
3 focuses on analysis of the source spectrum. Estimation of the intrinsic emission of the 
central engine (i.e., corrected for absorption along line of sight) is discussed in section 4. 
X-ray properties of IC 2560 and other extragalactic H2O maser sources believed to originate 
in accretion disks are considered in section 5; there we put IC 2560 into context and we 
discuss possible trends that may indicate connection between characteristics estimated from 
VLBI observations of masers and X-ray observations of their host AGN. Throughout the 
paper we adopt Hq=70 km s" 1 Mpc" 1 . 
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2. DATA REDUCTION AND ANALYSIS 

IC 2560 was observed using the XMM-Newton EPIC instrument on 2003 December 
2. Data reduction was carried out using SAS version 7.0.0. PN and MOS datasets were 
processed through the standard pipeline, using updated CCF files and SAS tasks epchain 
and emchain, respectively. Light curves extracted from the data do not indicate any of the 
cameras observed flares and hence no time intervals required filtering. The final exposure 
time was 72.2 ks for the PN camera, 80.5 ks for MOS1 and 80.6 ks for MOS2. 

We extracted the source spectrum from a circular patch, 15" (2.95 kpc) in radius, 
centered on the source coordinates in the binned image (binsize 80). Tests using the SAS task 
epatplot for PN data show good agreement between the expected and observed distribution 
of counts for the single, double, triple and quadruple event patterns. The source spectrum 
is thus unaffected by pileup. We used a nearby source-free circular region of radius 40" to 
extract the background spectrum, selected so that the source and background are equidistant 
from the readout node. After normalizing for area, the PN spectrum had a 3.1% contribution 
from the background, while the MOS1 and MOS2 spectra show a lower contribution from 
the background, approximately 2%. The source and background spectra were binned so that 
each bin contained at least 20 counts. Models were fitted to the spectrum between energies 
0.3-13 keV. 



The host galaxy has a semi- major axis of approximately 90" ( Ide Vacouluers et al.\\1991\ ). 
We used archival Chandra data to quantify the contamination of the XMM spectrum, due 
to large-scale emi ssion from the h ost g alaxy. The Chandra image (from the 2005 observa- 



tion discussed by iMadejski et al\ (120061 ) , Figure [T]) shows X-ray emission from two discrete 
sources, approximately 9.5" and 6.1" from the brightest pixel (within 0.8" of the optical cen- 
ter, the offset is comparable to the positional accuracy of Chandra at 99% limit and slightly 
higher than the 90% limit of 0.6"). The discrete emission appears to be dominated by soft 
and medium energy X-rays (< 2.5 keV). Diffuse emission in the plane of the galaxy is also 
visible. A circular region of radius 15" (size of the extraction region for the PN and MOS 
spectra), centered on the source yields 2201 counts, and 601 counts excluding the central 
4" for energies between 0.3 keV and 10 keV. The Chandra background as measured from a 
15" circle, placed approximately 90" from the center (comparable to the host galaxy semi- 
major axis) contains 435 counts. Thus, after background subtraction, we find the large scale 
emission from the host galaxy contributes approximately 9% to the source spectrum in the 
XMM observation. 
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SPECTRAL FITTING 



All spectral fitting was carried out using XSPEC version 11.3 flArnaud 1 11996). The 



spect rum shows all the key features observed in the Chandra spectrum (IMadejski et al 



20061 ) i.e. a hard continuum above 3 keV with a strong Fe Ka fluorescence line at 6.4 
keV, both indicative of a reflection origin in material far from the black hole. Below 2 keV 
there is additional soft emission, with multiple lines superimposed, that may be due to a 
complex mix of photo ionized and/or collisionally ionized plasma. We use the best-fit model of 
Madejski et all (120061 ) as a first approximation, consisting of two MeKaL collisionally ionized 
plasmas to describe the soft X-ray emissio n, a reflected component of th e (hidden) intrinsic 
power-law emission modeled by PEXRAV (IMagdziarz fe Zdziarski II1995T ) at higher energies, 
with a fixed inclination angle of 60°, together with the associated neutral fluorescence lines 
from Fe, Si and S Ka with rest frame energies of 6.4, 1.75 keV and 2.3 keV respectively 
(IBearden 1119671 ). We refer to this as Chandra-1, and use this as a starting point to fit the 
PN+MOS EPIC data from XMM- Newton (Tabled]). Fitting this to the data gives a good 
description to the overall spectral shape, but the better signal-to-noise XMM-Newton data 
show significant line-like residuals (Figure [2]). The strongest of these are at ~ 7 and 1.2 keV, 
corresponding to neutral Fe K(3 and Mg Ka, respectively, but there are also significant 
detections of neutral Ka emission from Ni and (marginally) Ca at 7.4 and 3.6 keV, as well 
as ionized Ka emission from He-like Fe at 6.7 keV. 

We also tested for the presence of other lines. Neutral Ka fluorescence lines from Ar 
(2.95 keV), O (0.55 keV) and Ne (0.84 keV) could also be present from the reflector. These 
are not significantly detected, but the up per limits t o the eq uivalent widths are consistent 
with expectations from solar abundances iMatt et all (Il997bl ). Similarly the data only give 
upper limits to the emission from other ionized lines which might be expected to accompany 
the He-like Fe at 6.7 keV such as H-like Fe Ka (6.955), He-like Ni Ka(7.8 keV) blended 
with He-like Fe K/9(7.88 keV) and H-like Ni Ka(8.09 keV). However, there is clear evidence 
for further ionized line emission at low energies, which is well matched by including a third 
plasma component. While this model (XMM-1, 3 MeKaL, 8 Gaussian lines and a PEXRAV 
continuum, Table [1]) matches the data fairly well (x 2 /d.o.f=358. 7/312), it is physically in- 
consistent. The He-like Fe line at 6.7 keV is either from a photoionized plasma, in which case 
it should be accompanied by a scattered power-law continuum, or from collisionally ionized 
plasma, in which case it should be produced by an additional MeKaL component. However, 
the required temperature of this 4th hot plasma emission is 4.4^'g keV (XMM-2, x 2 /d.o.f.= 
322.0/311, Table[T]), which seems high for starburst emission . It could indicate a violent 



interaction between the interstellar medium and an outflow (Veilleux fe Bland-Hawthorn 



19971 ). however, the detected outflow velocity (ISchulz fe Henkel 1120031 ) of 100 km s 1 is fairly 
low. Instead we favor the scattered power-law for a photoionized origin of the 6.7 keV line, 
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and note the resulting (marginally) better x 2 /d.o.f = 317.4/311. This gives our best-fit 
model XMM-bf (Figure ED, to the XMM PN+MOS data, as described in Table 1. 



3.1. Line Emission 



The lines that contribute significantly to the spectrum and their equivalent widths are 
summarized in Table [2j The equivalent widths for mostly neutral species are computed with 
respect to the reflected power-law compon ent. The spectrum h as no t evolved subst a ntiall y 
since the Chandra observations reported by llwasawa et all (120021 ) and lMadejski et all ( 120061 ). 
Specifically, we do not detect any statistically significant variability in the continuum flux 
in the 2-10 keV range or in the equivalent widths of the three emission lines detected in 
the Chandra dataset (Table [3]), though S/N is limited. The best-fit value for the 6.4 Fe K 
a line in the XMM dataset for model Chandra-1 is close to the lower limit for the same 
in the Chandra dataset. We attribute this to blending of 6.4 and 6.7 keV emission lines 
in the Chandra observation. The Fe Ka line in the XMM spectrum is marginally resolved 
and favors an intrinsic width of ~20 eV and <50 eV, similar to the Chandra data, and 
could originate in material at a few thousand km s _1 , typical of the BLR/inner torus. The 
emission features at 1.25 keV, 3.6 keV, 6.7 keV, 7.05 keV and 7.47 keV as well as the scattered 
continuum component are absent from the Chandra model, but these were not required to 
obtain a reasonable fit to the earlier spectrum. 

We inspected the spectrum for the presence of a 6.33 keV Compton down-scattered 
shoulder to the 6.4 keV line. To constrain this, we focus on the data above 5 keV, including 
only the lines within this energy range, together with the PEXRAV and scattered power 



law as continuum components with phot on index fixed 



Dut not norms) to the best-fit value 



given in Table 1. We included a m odel (jMadeiski et al 



Compton down-scattered shoulder (jlllarionov et a/.lll979l ) and fit this to the data but obtain 



20061 ) to describe the shape of the 



only an upper limit of about 10 per cent of the narrow Ka 6.4 keV line. This is at the l o w end 
of what is expected from reflection, indicating a face on vi ew of the reflector (IMatt I (120021 ) 
assuming the geometry described by lGhisellini et all (119941 )) and is in contradiction with the 
edge-on view of the disk inferred from maser emission. The face on view is also inconsistent 
with the expected orientation for a Seyfert 2 source under the unification paradigm. 

A similarly low Co mpton shoulder was observed for another Seyfert 2 source, Mrk 3, by 
Pounds fe Page I (120051 ) . However, they point out that this may be due to an error in the 
registered energy of split-pixel photons, which can be overestimated by up to 20 eV. Following 
their procedure and selecting only single-pixel events, we fit a Compton shoulder of 1.9l^ 2 9 5 % 
for a simultaneous fit to the PN and MOS data. A slightly higher limit was measured at 
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.7 



+16. 2c 
-8.9 ' 



for PN data alone, though this is probably due to decreased signal-to-noise. 



The fitted Compton shoulder height should be viewed with some caution. Detection of 
the shoulder is close to the limit of what is possible with available data. The inferred intrinsic 
width for the Fe Ka line is comparable to the shift in split pixel registration discussed above. 
The line width may also affect the detection of the Compton shoulder, due to limitations 
in modeling. Hence we adopt an upper limit of 25% on the Compton shoulder but stress 
that the data at face value are inconsistent with a reflection (or transmission) origin. On 
the other hand, recent work by lAwaki et all (120071 ) shows 10% upper limit for Mrk 3 from 
higher resolution Suzaku data, consistent with the XMM data at face value. IC 2560 may 
have a similarly lo w Compton shoulder, due to a clumpy torus structure as suggested by 



Awaki et all (120071 ) . However in that scenario, it will be very difficult to produce the observed 



amplitude of the reflected component. 



3.2. The Soft X-ray emission 



The soft X-ray emission dominant below 2 keV is almost certainly due to hot gas. For 
modeling purposes, we have used three MeKaL components to obtain a good description of 
this emission. This gives plasma temperatures of 0.1 keV, 0.32 keV and 0.66 keV. Replacing 
the three components with a single model where the plasma shows a power-law gradient with 
index 0.001 (best-fit value) in temperature does not give a satisfactory description of the low 
energy emission (% 2 /d.o.f = 616.8/315). However, some or even all of the observed emission 
below 2 keV could also be due to the photoionized scattering medium. Non-nuclear sources 
in the host galaxy could also have significant contribution t o the observed spectru m below 2 
keV. We note that the spectral analysis of Chandra data by iMadejski et all (120061 ) indicated 
that a collisionally ionized plasma was required even in presence of photoionized gas, in 
order attain a realistic description of the spectrum. However, the analysis neglected the 
scattered power-law continuum that must accompany any such photoionized emission, and 
the complexities of radiativ e transfer in the 6.7 keV line which can substantially change the 
strength of the line feature (IBianchi fe Matt 1120021 ) . A pure photoionized origin is suggested 
by the detection of radiative recombination line s associated with OVII at 0.69 keV in RGS 
data for this source (IGuainazzi fe Bianchi 1120071 ). though the signal-to- noise in these data is 
very limite d. The dominance of pho toionization in similar Compton th ick sources such as 
NGC 1068 (IKinkhabwala et q/.ll2002l ) and Circinus (IMassaro et a/.ll2006l ) implies it may be 
an important mechanism in IC 2560 as well. However, given the high density of line emission 
and low spectral resolution at these low energies, it is difficult to disentangle collisional and 
photoionization components and measure their relative contributions to the total emission. 
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ESTIMATE OF UNABSORBED EMISSION 



The obscurat i on co lumn density for IC 2560 is very high, > 10 24 cm -2 , as inferred 
by llwasawa et all (120021 ) and iMadejski et al\ (120061 ) . This so urce has not been dete cted at 
higher energie s (above 10 keV) by th e RXTE Slew Survey (jRevnivtsev et al\ 120041 ) or by 
SWIFT/BAT jMarkwardt et q/Jl2005h . If we extrapolate the best-fit spectrum to 200 keV, 
we find that the total flux from the continuum components in our model would contribute 
7.8 xlO -13 erg cm~ 2 s" 1 in the energy range 3-20 keV relevant to the RXTE Slew survey 
and 2.8 xlO -12 erg cm -2 s -1 in the energy range 14-195 keV corresponding to SWIFT/BAT. 
These values are not far below the detection limit of 10 -11 erg cm -2 s _1 for both surveys, 
strongly limiting the amount of direct continuum in these bands from the central engine. 



The [OIII] luminosity (L[q tjt]) can be used to g a uge the (unabsor bed) intrinsic 2-10 keV 
luminosity (Li2_ 10 ) of an AGN (IPanessa et a/.ll2006l ). iGu et al\ (120061 ) use stellar population 
synthesis modeling to estimate and subtract the stellar component from observed spectra, to 
predict a nuclear [OIII] flux of 1.25±0.05xl0 -13 erg cm -2 s _1 , after a redde ning correction 



of <1% determined from and Hg line ratios. We refit the data from IPanessa et al. 



( 120061 ). excluding the points with only upper limits, and scaled the parameter uncertainties 
by \/(x 2 /d.o.f) to obtain Logi V 2 _ 1Q = (1.16±0.04)(Log 10 L om -40.22) + (41.57±0.05). This 
gives us an approximate Yj\_ 
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6.17lg;|gxl0 41 erg s" 1 . 



In the following discussion, intrinsic emission refers to the flux or luminosity (L2-10) 
of the central engine that would be visible if there were no obscuration, direct component 
corresponds to intrinsic emission visible after it passes through the obscuring column suffering 
a reduction in low energy flux. We can also obtain an o rder of magnitude est imate for the 
intrinsic luminosity from the IR luminosity of the source, llwasawa et al\ (120021 ) estimate the 
Ljji from the IRAS to be 3xl0 43 erg s _1 . Assuming a power-law of photon index 1.8 from 
the disc at lOeV to lOOkeV and that all the observed IR emission is due to reprocessing 
of the X-ray emission, this gives intrinsic luminosity L2_io=6xlO erg s _1 . Stellar heating 
may have a significant contribution to the observed IR flux. The L2_ 10 from IR and OIII 
measurements are within an order of magnitude which is acceptable in light of uncertainties 
in both estimates. 

The estimated L2_ 10 is consistent with an unabsorbed power-law with photon index 
1.75 (best-fit value in XMM-bf) and normalization of 1.3xl0~ 3 . From this we estimate the 
reflected solid angle (f2/27r) from neutral material ~ 0.75 assuming the material is viewed at 
60° inclination, and the scattered fraction from the photoionized material to be ~ 0.006. We 
use th is intrinsic power la w to estimate the obscuring column density. For densities > 10 24 



cm 



Matt et al\ ( 11999a! ) show that electron absorption due to Compton scattering can be 



significant and photoelectric absorption alone is not sufficient. For a covering fraction <Cl, 
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electron scattering out of the line-of-sight is the only contributing process, and absorption 
can be sufficiently modeled by a combination of XSPEC models "wabs" and "cabs" . The 
unabsorbed power-law parameters were frozen to the values described above while the column 
density was varied and x 2 compared with the best-fit value. We found that at column 
densities less than 3xl0 24 cm -2 , the A% 2 due to the presence of the direct component is 
>2.78, so it contributes noticeably to the observed XMM spectrum. Thus the column must 
be > 3xl0 24 cm -2 . This is so modified by electron scattering that it makes little difference 
to the extrapolated SWIFT bandpass, so that the total continuum flux (direct + reflected + 
scattered components; FigureH]) in the 14-195 keV regime is 3x 1CT 12 erg cm -2 s" 1 . However, 
the reflector has a large solid angle as inferred from the PEXRAV modeling. In this scenario, 
electron scattering into the line-of-sight might partially compensate for scattering out of the 
line-of-sight and the required absorption column density could be higher. For a covering 



fract ion of 1, the direct component may be modeled using PLCABS, developed by (jYaqoob 
19971 ). For a column density of order 10 24 cm -2 , the inferred flux for this model, in the 



14-195 keV regime is fairly similar. The difference, and contribution of electron scattering 
to absorption, becomes more significant for higher column densities (Figure [5]). Substituting 
PLCABS to model the direct component with an absorbing column of 3xl0 24 cm -2 , gives 
a somewhat higher flux of 4xl0~ 12 erg cm~ 2 s^ 1 in the 14-195 keV regime. The flux values 
from both models are significantly lower than the SWIFT/BAT sensitivity, consistent with 
the non-detection of IC 2560. 



5. CONSTRAINING THE SOURCE GEOMETRY 

5.1. Warm Scatterer 

The detection of line emission at 6.7 keV due to He-like Fe (FeXXV) indicates presence of 
a warm, partly-ionized medium. The 6.7 keV feature is intrinsically a composite of 4 emission 
lines, unresolved in the XMM data: a forbidden line at 6. 65 keV, two interc ombination lines 



at 6.67 keV and 6.68 keV and a resonant line at 6.7 keV (IMatt et a/.lll996l ). EPIC does not 
have sufficient energy resolution to differentiate between these energies. Nevertheless, when 
allowed to vary, the rest frame energy of the line in the spectral model converged to 6.7io2 
keV. This may indicate that the resonance feature marginally dominates the to tal emission 
of the complex and implies densities less than 10 23 cm -2 , in accordance with IMatt et al. 



(119961 ). much less than that for the central source. We conclude the scatterer and central 
source are well enough separated that we have a relatively clear line-of-sight to the former. 
We associate this warm gas with the scattered power-law component that dominates the 
model continuum around 2 keV. In the energy range 2-10 keV, this component contributes 
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6.7^' 2 xl0 39 erg s _1 or ~1% of the the intrinsic emission incident on the scattering medium. 
We can estimate a lower limit to the column density of the scattering medium from the 
fraction of scattered flux, 

I-> scattered ^ 



-NhCTt 



T 



-^Intrinsic 4"7T 

where Q/Aw is the covering fraction, N# is the column density and ot is the Thomson cross- 
section. We find, (Q/Att) Nh"=1.5x10 22 cm -2 for the scattering medium. For a covering 
fraction < 1, this is a lower limit for the scattering medium. 



We can apply the analysis of iBianchi fc Matt I (120021 ) to obtain some order of magnitude 
estimates of the ionization parameter and distance of the warm scatterer from the central 
source. The ionization parameter is defined to be 



U x 



J2keV 



(2) 



Aimer 2 

where L is the ionizing luminosity, n the particle density and r the distance of the scattering 
cloud from the source. The equivalent width of the 6.7 keV line, 2.011};q7 keV measured 
with respect to the scattered continu um is significantly high er than expected for the above 



range in estimated column densit y (IBianchi &: 



turbulent motion of the warm gas (IBianchi et al 



Matt 



2002j), unless we take into account 



2005al ). The turbulent motion of the scat- 



tering material will result in line broadening, giving rise to larger obs erved equivalent width. 
Nonetheless, as can be seen from Fig. 4 in IBianchi fc Matt I (120021 ). the equivalent width 
is underestimated for all densities, hence adopting the ionization parameter corresponding 
the the peak in equivalent width for each column density, we estimate XJ X ~0.32. Though 
we do not detect the FeXXVI line at 6.955 keV, we obtain an upper limit of 196 eV from 
the current data. This yields a lower limit on of 0.2. For a photon index of 1.75, Equa- 
tion 2 reduces to U x = 0.23(L 2 _i /(W 2 )). We can express the particle density in terms 
of the column density of the scatterer and the cloud size R, nr 2 ~ t 2 (Nh/R) from which 
we obtain L 2 -io/£4 = r 2 (Nh / (0. 23 R)). If we assume the cloud size is comparable to its 
distance from the center, for a covering fraction of 1, r is on the order of lOpc. For a smaller 
covering fraction, the column density of the scattering cloud will be higher and distance of 
the scatterer to the central engine will be smaller. Nonetheless, our characterization of the 
warm scattering medium is subject to some important caveats. The ionization parameter 
is particularly uncertain. Specifically, we have assumed the relationship between equivalent 
width and ionization parameter will be relatively unchanged despite the introduction of a 
turbulent velocity, and rely on a scaling factor to account for the higher observed equivalent 
width. 



We note that the observed X-ray spectrum of IC 2560 does not show any intermediate 
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ionization lines, i.e. between the mostly neutral Ka lines due to cold Mg, Si, S, Ca, Fe 
and Ni, and the highly ionized He-like Fe at significantly higher temperatures. Emission 
lines from diverse ionization states of these and other elements have been detected in other 



Seyfe rt 2 sources such as Circinus (IMassaro et a/.ll2006l ) and NGC 1068 fjPounds fe Vaughn 



20061 ). It is likely that we lack the sensitivity to detect these lines. It is also possible that any 



intermediate ionization lines are blended with the detected Ka lines. This would provide 
a natural explanation for the relatively high observed equivalent widths. Source geometry 
might provide a third scenario where the transition region between the warm scatterer and 
cold reflector is either obscured from our line-of-sight or the scattering and reflecting media 
are (spatially) well-separated. An optically thick, mostly neutral medium is responsible for 
the reflection of intrinsic emission, observed above 3 keV while an optically thin, mostly 
ionized medium, with significantly lower column density scatters intrinsic emission along our 
line-of-sight. An obscured transition region between the two media, with no detected contri- 
bution to the observed spectrum would require a very special line-of-sight. Alternately, the 
scattering medium may reside above the accretion disk, as would a corona or an outflowing 
wind. 



5.2. Reflecting medium and Covering fraction 



The presence of Ka emission from mostly neutral species in the reflected emission in- 
dicates a low degree of ionization of the reflecting medium. We hypothesize that the cold 
reflecting medium and the obscuring material along our line-of-sight are one and the same. 
In that case, the cold reflecting medium has a very high column density and is optically 
thick. The fluorescence lines detected in the XMM spectrum probably originate in the re- 
flecting medium and aid in better understanding its composition and geom etry. The Fe Ka 
line width is consistent with sup ersolar abun dance (IBallantyne et aZ.I |2002| ) as is the upper 
limit for the Compton shoulder fjMatt II2002I ). The obser ved line widths for the remaining 
fluorescent lines are also much higher than predicted by iMatt et al\ (Il997al ) for a 60° and 
90° inclination angles, in most cases, by a factor of 2 or more. Th is may impl y a met al-rich 
environment a nd /or a different g eometry than that discussed by IMatt et all ( 11997a! ) . The 
calculations in IMatt et al\ (I1997al ) are most relevant for high angles of incidence and a plane 
parallel slab illuminated by an external source. A different geometry of source illumination 
may affect the observed equivalent widths and may explain the high values obtained for lines 
observed in IC 2560. We cannot rule out chemical enrichment as an alternative. The observed 
equivalent widths are consistent with 2xZ . However, supersolar abunda nces of lower met- 



als w ould result in partial absorption of the higher energy Fe Ka photons fjGeorge fe Fabian 
199ll ). This implies that the Fe abundance could be ~3-4xZ . Enhancement of low energy 
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Ka line widths due to blending with L-shell Fe lines is also likely. In this case, it is possible 
that the lower elements have solar metallicities and Fe abundance is ~2xZ . 

The 25 % upper limit on the Compton shoulder obtained from the comp osite PN an d 
MOS data is consistent with a range of inclination angles, from 30° to 90° ( iMatt 1 120021 ). 
We propose three different geometries for the reflector based on current data, where the 
reflector could be a warped thin disk, a thick clumpy disk or a flared disk inside a clumpy 
torus (Figure [6]). The distribution of maser e mission is consistent with an edge-on view of 
a thin Keplerian disk with little or no warp (jlwasawa et al\\2002\ ). If the maser emission 
originates in a relatively flat part of the disk, or our line-of-sight is such that we do not 
observe the warp, for e.g. orthogonal to a line of nodes, the warped disk scenario may be 
viable. However, it would be difficult to produce the relatively large (13%) fraction of the 
intrinsic emission observed in reflection, from a warped disk. Also the covering fraction for 
the obscuring medium inferred from modeling (0.75, refer section 4) would imply an extreme 
warp. Given the special constraints required to reconcile observations, a warped disk seems 
unlikely. A thick clumpy disk, w ith a filling factor high enough to obscure the central engine 
may be a feasible alternative. iNenkova et al\ (120021 ) show a clumpy torus structure can 
produce the necessary obscuration to shield a Seyfert 1 nucleus and produce a Compton 
thick source. In this case the masers would trace the equatorial plane of the disk toward the 
inner radi us, where pump energy from th e central engine can penetrate the clumpy disk as 
shown by iKartje. Konigl fe Elitzurl (119991 ) in the particular case of NGC 1068. The width 
of the Fe-line would be dictated by virial rotation. The observed 30 eV intrinsic width is 
consistent with the estimated central engine dynamical mass of 2.9 xlO 6 M and a radius of 
~ 0.04 pc, assuming Keplerian rotation. Radiative forces may drive a thi rd scenario with 
a thin inner disk and inflated outer disk (jChang et al\ 120071 ; iKrolik 1 120071 ) which can also 
produce the observed reflection fraction, equivalent widths and the intrinsic Fe Ka width. 
The later would be determined by the somewhat slower rotation of the torus and associated 
turbulent motions in the thickened disk, due to radiative heating, which can be comparable 
to gravitational forces at these radii. 



5.3. Comparison with other H 2 disk maser host galaxies 

We compare IC 2560 to the existing (though small in number) sample of known maser 
hosts with similar sub-parsec scale edge-on disk structure. The selection criteria were: (1) 
maser emission originating in a disk as determined from VLBI studies or inferred from maser 
spectra featuring red and blueshifted, high velocity Doppler components (so that the mass 
of the black hole can be estimated); (2) high inferred obscuring column density, in excess 
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of 10 24 cm 2 . We also include NGC 4258 as this is the archetypal maser host galaxy even 
though it has column of < 10 23 cm -2 in contrast to the Compton thick objects in our sample. 

Table 4 gives observed details of the eight sources in the sample. IC 2560 is typical of 
these in many respects such as equivalent width of the Fe K line (NGC 3079 and Circinus), 
reflected fraction (Mrk 3, Circinus) and Compton shoulder (NGC 1068, Mrk 3 and Circi- 
nus). All this suggests a similar geometry for the reflecting and absorbing material in these 
source s (edg e-on geometry where the obscuring/reflecting material is along line of sight; 



Matt I (120021 )). For a zeroth order analysis where we assume thermodynamic equilibrium for 
these sources, the temperature should be Ai:R 2 aT A ~ Lboi (both sides also depend in the 
same way on the covering factor of the material, Q/4tt, as this controls both the emission area 
and the fraction of illumination intercepted). Thus T 4 oc L^/ R 2 oc L 2 _ 10 /-R 2 for a constant 
bolometric correction, which allows us to crudely estimate the temperature for the region of 
the disk where the maser action occurs. For the observed Ll_ l0keV /L Edd taken from Table 
[5j we assume a bolometric corre ction factor of 20, consistent with the range 15-25 suggested 



by I Vasudevan fe Fabian I (120071 ). 



The significance of this relation to the data is limited by the large uncertainties in radius 
and luminosity. Disk geometry (e.g., warping, inclination to the line of sight) and maser 
beaming can substantively affect the inferred inner radius. We assume an uncertainty of 20% 
in the radius due to these factors and use it in weighting the rms scatter. The luminosity 
estimates depend on complicated model fitting to spectra, which can be particularly difficult 
for the high obscuration seen in these maser hosts. We cannot quantify the uncertainty in 
luminosity and instead note that the following analysis should be treated with due caution. 

A limited fit to the data in log-log space can be obtained by including uncertainties in 
radius alone and yields a curve with best-fit temperature ~600 K. Figure [7] shows R plotted 
against L hoh together with the predicted relation for constant temperatures of 300 K, 400 K, 
500 K, 600 K, 700 K, 800 K, 900 K and 1000 K. The observed inner radii of maser emission 
may be consistent with temperatures in the range 500-700 K, this is in ag reement with 



the predicted temperature range of 400-1000 K for maser disks (jElitzur 1119921 ). In physical 



disks, inefficiency of radiative heating is likely balanced by viscou s heating and rad iative 



dissipation, which have not been included in this simplistic analysis (IDesch et a/.lll998l ). The 



temperatures may be even higher for the more luminous sources whe re the Eddington ratio 



may be more consistent with a higher bolometric correction of 40-70 (jVasudevan fe Fabian 



20071 ). Conversely, if the estimated unabsorbed luminosities for the sources in this sample 
are systematically lower, they will imply smaller disk temperatures. Better understanding 
of luminosity estimates in this sample and the associated uncertainties are necessary for a 
more detailed analysis. 
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Warps are a marked feature of some of the accretion disks traced by water masers, and 
the degree of warp may also show some correlation with X-ray luminosity. We divide the 
sources into two groups, A and B, depending on the distribution of maser emission in the 
disk. Group A includes IC 2560, Circinus and NGC 3393 all of which show maser emission 
along a thin disk in the equatorial plane and Keplerian rotation curves, similar to NGC 
4258. Group B consists of NGC 1068, NGC 4945, NGC 3079 and NGC 1386 all of which 
show maser emission from a geometrically thick structure with flattened rotation curves. 
We note that the largest warp in Group A appears in the highest X-ray luminosity galaxy, 



Circinus ( IGreenhill et al.ll2003l ). IC 2560 and NGC 3393 show X-ray luminosities lower by at 
least a factor of 2 while we infer that the warps are smaller by factors of < 3, from maps of 
maser distributions in these two sources (Greenhill et ai, 2008, in prep). NGC 4258 exhibits 
variation in Euler angles that is several times smaller than Circinus for a X-ray luminosity 
lower by an order of magnitude (Table [5]). These observations may be consiste nt with models 



in which radiative to rques driven by disk irradiation are responsible for warps (jPringle 111996 



Maloney et a/.lll996l ). On the other hand, the absence of strong warps among AGN in Group 
B is difficult to understand, we speculate that presence of clumpy disk structure in Group 
B and the associated large variations in optical to infrared re-radiation may reduce torques 
that drive the orbits of disk material out of the plane and warp the disk. 

We charac terize the degree of substructure within the disk using the Toomre Q-parameter 



( iToomre II1964I ). It is also known as a disk stability criterion and can be used as an indicator 
of possibility of fragmentation due to self-gravity. For Q^>1, the impact of self-gravity is 
small and implies stability of a smooth distribution of disk material. For Q<1, the reverse is 
true, and the disk material may be anticipated to clump despite orbital motion in a central 
potential. Assuming uniform gas disks with observed central mass, at a given distance R, Q 
is defined as 

where fl is the angular velocity, v c is velocity dispersion and E is the surface mass density. We 
approximate the dispersion using the mean velocity width spanned by maser clumps along 
the rotation curves traced in position velocity diagrams. For a disk in pressure equilibrium, 
the surface density is <n>rriH 2 v c /v rot R, where v rot is orbital velocity at the radius R, 
<n> is the number density of the disk gas and is molecular mass. There is significant 
uncertainty (up to an order of magnitude) associated with the Q parameter discussed here, 
mainly due to the range of densities that can support maser action (10 8 -10 10 cm -3 ); we adopt 
10 9 cm -3 . We find that at the inner radius of maser emission, disks in Group A have a Q>1 
(Table [5]). A similar calculation of the Q-parameter is not appropriate for disks in Group 
B where maser observations provide significant evidence of self-gravity in the disk (for e.g., 
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NGC 3079: kondratko et all fl2005h . NGC 1068: lLodato fe Bertin I (J2003J)) 



If the assumptions we make in computing the Q parameter are broadly valid, the range 
spanned by Q-parameters for these sources may point to a u nified picture of accretion disk 



structure in maser hosts. Assuming hydrostatic equilibrium, iMilosavljevic fe Loeb I (120041 ) 
suggest a characteristic distance scale of 0.42pc for a central mass of the order of 10 7 M® 
beyond which the disk may be marginally unstable and likely to form clumpy structures due 
to self-gravity. In the unified picture, the overall disk structure in all these sources could 
be the same, with the maser emission in Group A sources tracing the inner, gravitationally 
stable parts of accretion disks and gravitationally unstable outer parts in Group B (Table 
Figure [7]). We focus on the inner radii of the maser disks, which we may a nticipate to be 



l inked to intrinsic luminosities via the temperature of dust sublimation, ~10 3 K f lNeufeld et al 



( 119941 ) : the outer radii of the maser disks, may be determined by the contribution of other 
factors such as density profiles, shielding and re-emission from surrounding gas). The lower 
luminosity of sources in Group A allows survival of molecules (and maser action) at smaller 
radii, where the disk structure is well-ordered. Among higher luminosity sources in Group B, 
molecular emission arises farther from the black hole, which is also where the disk structure 
is likely to be clumpy. X-ray observations of other Compton thick maser hosts, where VLBI 
data points to a disk origin will help expand the sample size and test the preliminary trend 
reported here. 



6. SUMMARY 

A deep 70 ks XMM-Newton spectrum for IC 2560 shows evidence for a warm scattered 
continuum component as well as number of emission lines not detected in earlier Chandra 
observations. These include Ka emission from Mg (1.25 keV), Ca (3.6 keV) and Ni (7.47 
keV) as well as Fe K/3 (7.05 keV) and Ka line emission from Fe XXV (6.7 keV). We do 
not detect line emission from species at intermediate ionization states, between the nearly 
neutral material that gives rise to the Ka emission, and the warm scattering medium with 
FeXXV. The direct X-ray emission from the central source appears to be completely obscured 
up to 10 keV; the inferred column density is in excess of a few times 10 24 cm -2 . This is in 
agreement with a 25% upper limit on the magnitude of the Compton shoulder. The reflecting 
medium is mostly neutral, optically thick and could also provide the obscuring column toward 
the central engine. The 6.7 keV line is associated with the optically thin, mostly ionized, 
scattering medium, with significantly lower density. The dominance of the resonance line at 
6.7 keV indicates the density has to be lower than 10 23 cm -2 , while the estimated fraction 
of intrinsic emission scattered along line-of-sight indicates a lower limit of 10 22 cm -2 . The 
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scattering medium can plausibly exist on the scale of a few parsec, but the exact geometry 
is unknown. Limits on density of the scatterer indicate a covering fraction between 0.1 
and 1. Using a small sample of H 2 disk maser hosts, selected using physically motivated 
criteria, we find that the unobscured luminosity of the central source may shape molecular 
disk structure in this sample. The inner radius as well as degree of substructure in the disk 
may be related to the unobscured X-ray emission, though a larger sample is required to 
confirm this trend. 

We thank E.M.L. Humphreys for many useful discussions and Julian Krolik for help- 
ful comments. This research has made extensive use of the NASA/IPAC Extragalactic 
Database(NED) which is operated by the Jet Propulsion Laboratory (JPL), California Insti- 
tute of Technology, under contract with NASA, and the NASA's Astrophysics Data System 
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Table 1. Characteristics of XMM EPIC PN+MOS spectrum 



Model name 


Model 


i _1_U CIIl 1 


kT[l] 

V * / 


kT[2] 

IKc V f 


kT[3] 

^Ke v J 


kT[4] 

V V J 


r 


Eq. W.( 6 > 

V * / 


X 2 /dof 


Chandra- 1 


ab(pexrav+2mekal+3gauss)M 




n fi o+0.01 


o.o8iS:gjw 






Q 1 O + 0.01 


2 4n+ 017 


493.2/317 


XMM-1 


ab(pcxrav+3mckal+8gauss)' e ' 


1 6 +0 ' 2 


n 77 +0.05 


o.osig:SJw 


n on + 0.03 
u ' oz -0.03 




Z - Z9 -0.13 


2 7n+ - 18 

z - ' u -0.18 


358.7/312 


XMM-2 


ab(pexrav+4mckal+8gauss) (f) 


14 7+0-2 


71 +0.06 


o.o8tg:SJ« 


n si+ 04 

u ' al -0.05 


44 +2.9 


1 62+ 021 


2 7n+ - 18 

z> ' u -0.18 


322.0/311 


XMM-bf 


ab(pexrav+powerlaw+3mekal+8ga 


,uss)(f) U.7±ll 


n 7 o+0.06 
u - <°-0.09 


0.08iS:gJ« 


01 +0.03 
u " r>1 -0.06 




i 7^+0.21 
x - ' d -0.23 


9 oo+O.lS 
z "-" z -0.17 


317.4/311 



( a ) The lower limit on the absorption column density fixed to Galactic absorption. 

( b ) Equivalent Width for Fe Ka line, with respect to the reflected continuum. 

(^Includes contributions from line emission corresponding to Si, S, and Fe Ka at energies 1.75 keV, 2.31 keV and 6.4 keV as well as continuum emission from 
the reflected component and 2 MeKaL plasmas. 

W Parameter value converged to lower limit. 

( e 'lncludes contributions from line emission corresponding to Mg, Si, S, Ca, Fe Ka, Fe K{3 and He-like Fe Ka at energies 1.25 keV, 1.75 keV, 2.31 keV, 3.6 
keV, 6.40 keV, 7.05 keV and 6.7 keV respectively as well as continuum emission from the reflected component and 3 MeKaL plasmas. 

(■^Includes contributions from line emission corresponding to Mg, Si, S, Ca, Fe and Ni Ka at energies 1.25 keV, 1.75 keV, 2.31 keV, 3.6 keV, 6.4 keV and 
7.47 keV; Fe K /3 at 7.05 keV as well as continuum emission from the reflected component and 4 MeKaL plasmas. 

(^Includes contributions from line emission corresponding to Mg, Si, S, Fe and Ni Ka at energies 1.25 keV, 1.75 keV, 2.31 keV, 6.4 keV, 7.47 keV; Fe K /3 at 
7.05 keV and He-like Fe Ka at 6.7 keV as well as continuum emission from the reflected component, 3 MeKaL plasmas and a scattered powerlaw, with photon 
index same as the reflection model. 



Table 2. Line Equivalent Widths for XMM and Chandra data 



Instrument Mg K a 

1.25 kcV< a > 
(keV) 


Si K a SKa 
1.75 kcV< a ' 2.31 keV< a > 
(keV) (keV) 


CsKa 
3.6 keVM 
(keV) 


Fe K a 
6.4 kcV< a > 
(keV) 


Fe K p NiKa 
7.05 keV< a > 7.47 keV< a ) 
(keV) (keV) 


He-like Fe K a 
6.7 keVW 
(keV) 


XMM PN+MOS< 6 > 0.55^0 24 


n 55+ 014 

u ' oo -0.15 


n q 4 +0.15 


n ns+ - 08 
U.U8_ 08 


9 o 9 +0.18 
z ' oz -0.16 


n 3n+ ' 02 


O.+0.02 
u -°- L -0.02 


2 ni +107 ( c ) 

z - ui -1.07 


Chandra ACIS-SW 


33+ 17 


34+0-W 
u -°*_0.20 




„ 77 +0.49 
z - ' ' -0.49 









' a )Rest frame energies. 

''''Equivalent widths for XMM data using best-fit model XMM-bf, described in Table 1, calculated with respect to reflected 
component. 

' c ' Equivalent width calculated with respect to scattered component of the continuum. 

( d ) Equivalent widths for Chandra data using best- fit model from iMadei ski et al. calculated with respect to reflected 

component. 
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Table 3. Emission Characteristics from Chandra and XMM-Newton Studies' - 1 



Energy range Emission Chandra XMM-Newton 



6.4 keV 


Equivalent Width 


2.77±°;£ keV 


2.32±g;^ keV 


6.25-6.45 keV 


Observed Flux 


1.35± ;g^xl0- 13 erg cm" 2 s" 1 


1.23lo;ogXlO- 13 erg cm" 2 s" 1 




Observed Luminosity 


2.85l?;?|xl0 40 ergs" 1 


2.68± J ;°Jxl0 40 ergs" 1 


2-10 keV 


Observed Flux 


3 - 84 io:46 xl0_13 er § cm ~ 2 s_1 


3.881° " xlO" 13 erg cm" 2 s" 1 




Observed Luminosity 


8 - 06 ^3:o5 xl ° 40 er § s_1 


8.15i ,; 31 xl0 40 erg s" 1 



("'Values and ranges from model Chandra-1, fitted to Chandra data, described bv lMadeiski et "all i2006l ). and 
XMM PN+MOS data as described in Table 1. 



Table 4. Observed Properties of Compton-thick Disk Maser Nuclei 



Source 


(a) 
(M ) 


(b) 
Disk size 
(pc) 


(c) 
(Mpc) 


(d) 
(cm -2 ) 


(e) 

L Obs (2-10) 
(erg s" 1 ) 


(f) 
Eq. W. 
(keV) 


(g) 

CS 
(96) 


References 


IC 2560 


2.9±0.6xl0 6 


0.08-0.27 


41.4 


> 3xl0 24 


8.2 xlO 40 


9 ni+0.18 
-0.16 


<25 (h > 


(a),(b)- 1. 

(Hi fol ft 1 ( cr 1— f~ ri it; nsi npr 


NGC 4945 


~1.4xl0 6 


~0.3 


3.7 


~4xl0 24 


3xl0 38 






(a),(b)-2. 

(d) -3. 

(e) -4. 

(f) -5. 


NGC 1068 


~2xl0 7 


0.65-1.1 


14.4 


>10 26 


5.2 xlO 40 


1.2 


-20 


(a.),(D)-b. 

(d)-7. 
( e ) (g-)-8 

(f)-9. 


Mrk 3 


~ 4.5x10 s 




57.7 


1.36i»; c «xl0 24 


1.4xl0 42 W 


fil+ 03 
u ' Di -0.05 


iol 9 6 


(a)-10. 
(d)-ll. 
(e),(f),(g)-12. 


Circinus 


1.7±0.3xl0 6 


0.11-0.4 


1.1 


4.3± '; 4 xl0 24 


~10 4 °0) 


9 9 tr+0.26 


20±3 


(a),(b)-13. 

(d) -14. 

(e) -15. 

(f) -16. 

(g) -17. 


NGC 3393 


3.1±0.2xl0 7 


0.17±0.02 


53.3 


~4xl0 24 


3.1^ xlO 40 


1.4±0.8 




(a),(b)-18 
(d),(e),(f)-19. 


NGC 3079 


2xl0 6 


~0.4 


17.3 


~10 25 


~10 4(J 


2 4+ 2 - 9 




(a),(b)-20. 
(d),(e),(f)-21. 


NGC 1386 


1.2xl0 6 


0.44-0.94 


17« 


>2.2xl0 24 


4.9±0.9xl0 39 


1 s+ - 4 




(a),(b)-22 
(d),(e),(f)-19. 



(a),(b)-23. 

NGC 4258< ! > 3.9±0.3xl0 7 0.12-0.28( m > 7.2 6.0-15 xlO 22 2.4-10.4xl0 39 <0.049 (d),(e),(f)-24. 



( c ) Luminosity distance used in calculation of mass, radius and luminosity for each source 

(a),C>)>(e)Uncertainties are quoted where they are known. With the exception of NGC 4258, the mass, radius and luminosity in each 



object are uncertain by uncertainties in distance, as well as any modeling uncertainties. 
( e ) Observed luminosity in 2-10 keV energy range. 

" )Eq. W. of Fe Ka line measured with respect to the reflected component of the spectrum. 

( 9 )Compton shoulder strength relative to the Fe Ka line. 

W Upper limit obtained by considering only Pattern events. 

W Calculated from quoted flux in lBianchi et al. | ll2005bh . 

(j)\Vilson, private communication. 

(*' Distance adopted from iHenkel et a! I d2005h 

("NGC 4258 is a Compton-thin source and has been included here only for comparison. 

Source luminosity range observed over 10 years, variation due to variability in absorption column density. 
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Table 5. Inferred Properties of Compton-thick Disk Maser Nuclei 



Source 


^2-10 


L 2-1()/ L -Edd 


R/R G (i,) 


Qi„ (c) 


Q out (d) 


Group ( e 




(erg s -1 ) 












NGC 4258 


4.2-17.4xl0 4 °( / > 


1-4 xlO' 5 


3.2-7.2X10 4 


144 


11 


A 


NGC 3393 


4xl0 41 ^ 


lxl0~ 4 


5.7±0.4xl0 4 




32 


A 


IC 2560 


e.sxio 41 ^' 


2xKT 3 


2.8-9.5X10 5 


36 


0.9 


A 


Circinus 


1.2xl0 42 « 


5xl0" 3 


6.7-25X10 5 


8 


0.2 


A 


NGC 1068 


2.3xl0 43 < J ) 


lxlO" 2 


3.4-5. 8xl0 5 






B 


NGC 4945 


3xl0 42(fe) 


2xKT 2 


2.2xl0 6 






B 


NGC 3079 


6-12xl0 42 W 


2-4xl0" 2 


2.1xl0 6 






B 


NGC 1386 


9.6-38. 2xl0 42 ( m 


6-24xl0~ 2 


3.7-7.9X10 6 






B 



'"'intrinsic (unabsorbed) luminosity in the 2-10 keV range, after accounting for local 
and Galactic absorption. 

'^Radius for maser emission, in gravitational units. 

(c),(d)Q parameter computed for the inner and outer radii of maser emission for uni- 
form gas disks assuming, particle density of 10 9 cm' 3 characteristic of maser action 
from H2O at 22 GHz, and velocity dispersion, v c =1.5xl0 3 km s _1 corresponding to the 
sound speed of molecular gas. 

( e ' Grouping based on disk morphology, Group A: Thin, fairly flat disks, Group B: 
Thick disks. 

( g ) Luminosity from reflection dominated modeling of Beppo-Sax data bv lSalvati et all 
lll997f), preferred over tr ansmission dominated modeling, based on analysis of PN data 
bv lMaiolino et all lll998l) . 

W Intrinsic luminosity quoted by dKondratko et alll2006al) , using PIMMS, does not 
take into account obscuration due to Compton scattering and may be significantly un- 
derestimated. 



References. — (f ) jFruscione et all j2005)). fg)-|Sarvati ef"all ^99^"). (hVthis paper , 
(i)- Smith fc Wilson I COOll l. ( p-IPounds fe Vaughn I fcOOrt ), (kl jcuainazzi et all 1 I2OOOI) , 
Ql- llvomoto et all [120011 1, fml jKondratko et all l l2006al) 



Fig. 1. — Chandra image of IC 2560. North is up. Energies 0.3-1.5keV represented by 
red, energies 1.5-2.5 keV by green and 2.5-8 keV energy represented by blue colors. Diffuse 
emission from galactic plane of IC 2560 is visible as are two discrete sources of emission, 
approximately 9.5" from the central source, about 50° from the north and 6.1" from the 
center, about -20° from the north. The image has been smoothed using a Gaussian function 
with radius 2 pixels, to suppress statistical noise. Coordinates indicate brightest pixel in 
image 
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Fig. 2. — Composite PN and MOS data, binned with 20 cts/bin, fit to the Chandra best-fit 
model (top) and A% 2 residual (bottom). The model is described in Table 1 as Chandra- 1. 
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Fig. 3. — Best-fit unfolded model and spectrum for XMM PN and MOS data, binned with 
20 cts/bin, summarized in Table 1, XMM-bf. Both PN and MOS unfolded data are shown 
by red crosses. 
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Fig. 4. — Continuum spectrum extrapolated to 200 keV, with contribution from ionized 
plasma, scattered and reflected components as well as absorbed intrinsic emission of source. 
Absorption column density for intrinsic emission set to 3xl0 24 cm -2 . 
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5 10 20 50 10C 

Energy 

Fig. 5. — Direct component of intrinsic emission through obscuration columns of density 
10 24 cm -2 (dashed lines) and 10 25 cm -2 (solid lines), with only photoelectric absorption 
(blue), with photoelectric as well as electron scattering out of line-of-sight but not back into 
line-of-sight, corresponding to a very low covering fraction(green), and combination of both 
absorption processes including scattering into line-of-sight for unit covering fraction (red). 
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Fie. 6. — Cartoons of possible geometries for AGN in relation to prospective H 2 maser 
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Fig. 7. — Correlation between observed inner radius of maser emission and modeled intrinsic 
bolometric luminosity. Data overplotted with the relation anticipated for external radiative 
disk heating and thermodynamic equilibrium, R oc L 5 for disk temperatures 300 K, 400 K, 
500 K, 600 K, 700 K, 800 K, 900 K and 1000 K. Disks in Group B dominate the higher end 
of the luminosity distribution. 



